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We study the effect of disorder on the London penetration depth in iron-based superconductors. 

The theory is based on a two-band model with quasi-two-dimensional Fermi surfaces, which allows 
for the coexistence region in the phase diagram between magnetic and superconducting states in 
the presence of intraband and interband scattering. Within the quasiclassical approximation we 
derive and solve Eilenberger’s equations, which include a weak external magnetic field, and provide 
analytical expressions for the penetration depth in the various limiting cases. A complete numerical 
analysis of the doping and temperature dependence of the London penetration depth reveals the 
crucial effect of disorder scattering, which is especially pronounced in the coexistence phase. The 
experimental implications of our results are discussed. 

PACS numbers: 74.70.Xa, 74.62.En, 74.25.N-, 74.25.Dw 


I. INTRODUCTION 

Measurements of the magnetic penetration depth Xp 
as a function of temperature, doping, magnetic field, and 
crystal orientation provide invaluable information about 
the nature of superconductivity and the symmetry of the 
underlying order parameter (see e.g. reviews jl!2j ). In 
a single-component clean s-wave BCS superconductor, 
with order parameter A and a fully gapped Fermi surface, 
the low-temperature behavior of the London penetration 
depth SXl(T) = Al(T) — A^(0) shows an exponential 
decrease with temperature 

SXl(T) 

Al(0) 



(FeSCs) is of special interest. These materials have multi¬ 
ple Fermi pockets with electron-like and hole-like disper¬ 
sion of carriers. Because of a delicate interplay between 
interactions in various pairing channels, superconductiv¬ 
ity in FeSCs emerges in close proximity to a spin-density- 
wave (SDW) order, and the superconducting (SC) crit¬ 
ical temperature T c has a dome-shaped dependence on 
doping, with the T c maximum near the onset of SDW or¬ 
der.® 11 It has been proposed^ that superconductivity in 
FeSCs is unconventional, with the order parameter hav¬ 
ing opposite signs on different Fermi sheets, and named 
s 16 symmetry. The latter emerges because SDW fluctua¬ 
tions increase interpocket interaction, which is attractive 
for s 16 gap symmetry, to a level where it overcomes in¬ 
trapocket repulsion. Likewise, SC fluctuations tend to 
increase the tendency towards SDW. 


Even though nonmagnetic disorder does not directly af¬ 
fect A, it does modify the value of Xl, which becomes 
X~f 2 {T) oc A(T)ertanh[A(T)/2T], where a is the normal- 
state conduct ivityP® 

In contrast to the s-wave case, d -wave symmetry of the 
order parameter with nodes on the Fermi surface trans¬ 
lates to power-law temperature dependence for the pen¬ 
etration deptlP 


6X l (T) T 
A l (0) “a’ 


( 2 ) 


The power exponent of the low-temperature behavior 
is very sensitive to disorder scattering, such that 5Xl 
crosses over to quadratic behavior, <5 Al(T)/Al(0) ~ 
(T/ A) 2 , below a certain temperature scale T*, which is 
determined by the concentration of strong scatterersP 7 
The dependence of the penetration depth on various 
parameters in the case of iron-pnictide superconductors 


The emergent complexity of FeSCs with competing 
superconducting and magnetic instabilities, which may 
coexist in a certain region of the phase diagram] ISES] 
leads to peculiar dependencies of the penetration depth. 
Early experiments in 122-materials, Co- and K-doped 
BaFe 2 As 2 , revealed that down to the lowest tempera¬ 
tures and in a wide range of dopings the T-dependence 
of 5Xl can be systematically fitted by SXl oc T 2 ES HT ZI 
contrast, in 1111-compounds such as SmFeAsOi_ x F x — 
and PrFeAsC 1 ^the penetration depth has an exponen¬ 
tial temperature dependence consistent with a gap with¬ 
out nodes; no appreciable effect of scattering was ob¬ 
served. At the same time, data on another 1111-material, 
LaFePO 20 pointed out that SXl(T) varies approximately 
linearly with T, strongly suggesting the presence of gap 
nodes in this compound. Since these initial reports, the 
London penetration depth has been measured systemat¬ 
ically in a variety of families of iron-pnictides and iron- 
chalcogenidesPSiSS Perhaps the most striking recent ob- 
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servation is a disorder-induced topological change of the 
superconducting gap structure, as revealed from the low- 
T behavior of SXl in BaFe 2 (Asi_ x P x ) 2 p^ Nonmagnetic 
defects were controllably introduced by electron irradia¬ 
tion, and it was found that the nodal state of P-doped 
BaFe 2 As 2 changes to a nodeless state with increasing dis¬ 
order. Moreover, under further irradiation, the gapped 
state evolves into a different gapless state, thus providing 
evidence of unconventional sign-changing s-wave super¬ 
conductivity. Such unusual sensitivity of the supercon¬ 
ducting gap structure to disorder scattering is a unique 
characteristic feature of FeSCs. 

Theoretical studies of the penetration depth in FeSCs 
were discussed in Refs. |33H37] for clean samples based 
on the band model. The effects of disorder on the phase 
diagram, including pair-breaking scattering, and on the 
penetration depth were investigated in Refs. RSITO] . We 
study the effect of disorder on Ai in a systematic way 
and analyze its behavior in the part of the phase diagram 
where the SC and SDW phases coexist. On the techni¬ 
cal side, we develop a formalism that enables us to study 
the doping and temperature evolution of the penetration 
depth in the whole parameter space of the phase diagram. 
Recently electron irradiation was used to intr oduce disor¬ 
der into FeSC systems in a controlled way! 3 ° |32 ^ Thus, 
our theory is relevant for the interpretation of existing 
and future experiments along this exciting direction. 

This paper is organized as follows: In Sec. II we present 
our model, discuss underlying approximations and as¬ 
sumptions, and we analyze the phase diagram of the 
FeSC compounds. In Sec. Ill we derive and solve qua- 
siclassical Eilenberger equations with emphasis on the 
coexistence of SC and SDW orders. We then apply that 
formalism to study the London penetration depth across 
the whole range of the phase diagram, and at different 
temperatures. In Sec. IV we summarize our findings and 
place our work in the context of future developments. 


II. MODEL AND APPROXIMATIONS 


In this section we introduce the minimal model for 
iron-based superconductors in which doping acts as a 
source of disorder and produces a region of coexistence 
between superconductivity and magnetism. Further¬ 
more, right from the outset, we consider the case of 
nonzero external magnetic field that acts on orbital elec¬ 
tron motion, but assumed to be weak enough not to 
affect spin. We discuss the ground state properties of 
this model in zero field within the quasiclassical approx¬ 
imation, which we use later to compute the penetration 
depth across a wide doping range from the coexistence 
region to the purely superconducting state. 


A. Model 

Following the discussion in Refs. mm, we consider 
a model with two cylindrical Fermi surfaces. One Fermi 
surface has electron-type and another one has hole-type 
excitations. We introduce the following eight-component 
spinor 


^(r) = ($(r), V'c(r), V>}(r), ^/(r)), (3) 

where ^ Q (r) = (^(r), V’<4( r )) (a = c,/) is a Gor’kov- 
Nambu spinor, and tp\ a (r) are the creation operators for 
the electron ( a=f ) and hole (a=c) fermionic excitations 
at point r in real space with a spin component a =ti- 
The full Hamiltonian for the problem at hand 

r o:/3 

consists of kinetic part and interactions 

H(r) = H 0 (r) + H mf (r). (5) 

In the limit of weak magnetic field the noninteracting 
Hamiltonian matrix [Ho( r )]a/3 can be compactly written 
as 


A a Z6 

Hq(y) = —£f 3 p3<j 0 H-A • Vr 3 p 0 a 0 . (6) 

me 

Here, fj, pi, and cq with i = 0,1,2,3 are sets of Pauli 
matrices acting correspondingly in the band, Gor’kov- 
Nambu, and spin spaces; fo,/5o,<7o are unit matrices, A(r) 
is the vector potential, £ = — V 2 /2m — p, and p is the 
chemical potential. 

Interactions between the quasiparticles on the electron- 
and hole-like Fermi surfaces lead to the development of 
superconducting and spin-density-wave orders. Within 
the mean-field theory approximation, the corresponding 
expression for the interaction part of the model Hamil¬ 
tonian H m f reads 

H mf = H a + H m , (7) 

H a = -Af 3 p2d 2 , Hm = ti/ 3 3 M • <x. (8) 


Here, A is the superconducting order parameter, while 
M is the spin-density-wave order parameter. We ex¬ 
plicitly assume that the magnetic field is weak enough 
so that we can ignore the spatial dependence of both 
SC and SDW fields [see discussion after Eq. (14) be¬ 
low]. Also note that A and M must be computed 
self-consistently; we will derive the corresponding equa¬ 
tions in what follows. We emphasize that within the 
model under consideration, we study the case of s 31 pair¬ 
ing, i.e. the superconducting order parameters on the 
electron-like and hole-like Fermi surfaces have opposite 
signs, A( c ) = —A^ = A. We also ignore the possible 
mismatch due to differences in the band occupations and 
effective masses between the two Fermi surfaces. 
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Let us now introduce a disorder potential. In what 
follows we consider two types of disorder scattering: the 
first type is intraband disorder with potential Go, which 
scatters quasiparticles within the same band, while the 
second type with potential U v accounts for interband 
scattering. Thus, in the basis ^ for disorder potential 
we write 


well since both supercond ucting an d magnetic correlation 
lengths greatly exceed - \ F ! 41 I 42 I 44 I 45 I 

The central object in the quasiclassical approach is the 
Eilenberger function'll 

£ w (R,n) = — / ^-T 3 p 3 a 0 ■ G(iw„,R,p). (12) 

7 TV J Z7T 


U(r) = ^2 [U 0 T 0 P 3 V 0 + UnTifo&o] <5(r - Ri), (9) 

i 


where the summation goes over impurity sites. We as¬ 
sume that concentration of impurities is a;i mp . 

We will treat the effects of disorder within the self- 
consistent Born approximation. Specifically, we intro¬ 
duce a single-particle Green’s function in the Matsub- 
ara representatiorP as a solution of the following matrix 
equations 


iu n ~H( ri)-E u (R) G(iu n , ri, r 2 ) = /, 
-iu> n -H( r 2 )-E u (R) G(iuj n ,r 1 ,r 2 ) = I, 


( 10 ) 


where oj n = irT(2n + l)fo/5o<To, T is a temperature, 
I = TopoaoS(r 1 — r 2 ), and R = (ri + r 2 )/2 is the center- 
of-mass coordinate. In order to write down an explicit 
expression for the self-energy E(iw n ,R), in addition to 
the center-of-mass coordinate we introduce the relative 
coordinate r = ri — r 2 , and consider the matrix Green’s 


function (10) as a function of R and r. Furthermore, we 


perform the Fourier transformation with respect to the 
relative coordinate r, and in what follows we consider 
the function G(iw n ,R, p). Then, assuming that disor¬ 
der is uncorrelated, upon averaging over various disorder 
configuration^ we find the following expression for the 
self-energy (hereafter h = c= 1): 


(R-) = — 

7 TV 


4T^ 

7 TV 


d 2 p 
(27r) 


d 2 p 

W? 


T 0 p 3 a 0 G(iuj ni R, p)fop 3 CTo 


2 T1P3CT0 G(iuj n , R, p)ti/ 5 3 (J 0 , 


( 11 ) 


where the cross terms oc UoU n vanish; v is the sin¬ 
gle particle density of states, T 0 = 7w:ri m p|{/o| 2 /4, and 
= 7rz^rci mp | LGr | 2 /4- Clearly, the fully self-consistent 
computation of the order parameters A and M along 
with the self-energy E U (R) is a challenging problem. 
However, this problem can be solved efficiently using the 
quasiclassical approach. 


B. Quasiclassical approximation 


To derive an equation for the Eilenberger function 
£7 W (R, n) in a weak external magnetic field, one needs to 
eliminate the single-particle dispersion via a series of al¬ 
gebraic manipulations (see Appendix A for details). Tak¬ 
ing into account that the relevant values of the quasipar¬ 
ticle momentum p are close to the Fermi momentum p F , 
so that p/m ss vf n, we find the following equation for 

Gui- 


«w„f 3 /3 3 <7o,£ w (R,n) - -ff m ff 3 p 3 <T 0 ,&,(R,n) 
E w (R)f 3 /5 3 <X(h £ W (R, n) 
ev F n ■ A(R)fo/3 3 <j 0 , ^ W (R, n 


+ 


+v F n ■ (-iVr)^(R, n) = 0, 


(13) 


where the square brackets denote a commutator. Since 
we consider the limit of a weak magnetic field, we can 
look for the solution of this equation by perturbation 
theory, namely 

g w (R,n)=^ 0) +^ 1) (R,n), (14) 


restricting ourselves to corrections linear in powers of the 
vector potential A. This is why we could neglect the 
dependence of A and M on R, since the corrections that 
render both order parameters spatially inhomogeneous 
are of the order of 0(A 2 )PNext, we discuss the solution 
of the Eilenberger equation in the spatially homogeneous 
case. 


C. Phase diagram 


In this section we first review the ground-state prop¬ 
erties of the model (5l) with disorder by taking the limit 
of A = 0 in Eq. (113 1 and considering a uniform system, 
VGu' 1 = 0. We thus have 


WnT 3 p 3 <7 0 


gL 0) 


(H m { + E w ) • f 3 p 3 d 0 , G £ 0) 


= 0. 
(15) 

Without loss of generality we choose the SDW magne¬ 
tization to be along the z-axis, M = Me z . Then, the 
solution of © has the following form: 

G i 0) = gu>T 3 P3&o - ifuToPid 2 - is u r 2 Po^3, (16) 


The quasiclassical approximation is justified when the 
characteristic quantities for the problem at hand vary 
significantly on length scales that are much longer than 
the Fermi wavelength X F . In the context of iron-based 
superconductors, the quasiclassical approximation works 


where the functions g^, fun and s u are determined by the 
solution of the following system of algebraic equations: 

*A g u = fui (uj n + ^Y^guf) , , . 

iMg u = Su, (w n + 2T t g u ). 
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FIG. 1: (Color online) Phase diagrams in the (ro,T) plane 
computed by solving the self-consistency equations ( |18| ) and 
the Eilenberger equations (171 for various ratios between in¬ 
traband and interband scattering rates. T c o denotes the su¬ 
perconducting critical temperature in a clean system. All 
plots are obtained assuming T a o = 327,0, and T a o is the criti¬ 
cal temperature for the SDW state. 


Here we introduced the total scattering rate T t = rb+rV. 
In addition, the functions g^, f u , and satisfy the nor¬ 
malization condition gt — s ti = 1- Subsequently, su¬ 
perconducting and SDW order parameters can be found 
from 


iM 

9m 


A 

7tT ^ Su>, 

uj n >0 


iA 

9sc 


A 

kt /«> 

LJ n >0 


(18) 


FIG. 2: (Color online) Width of the coexistence region 

(r^ sdw) — Fg SC ^)/27rT c o is shown as function of F^/ro. The data 
points are found from the solution of Eqns. (181 and (171 at 
zero temperature. The plot is obtained assuming T a o = 377o- 



where g ac and g m are the coupling constants, and A is an 
ultraviolet cutoff. In the clean system, there is a phase 
transition from the paramagnetic to the SDW state at 
critical temperature T s o = 1.13Ae -2 /" 9m provided that 
gm > 9sc- If 9sc > 9m-, the ground state is a supercon¬ 
ductor with a critical temperature T c o = \.13Ke~ 2 ^ 9ac . 
We consider T s0 > T c0 , so that without disorder, the 
SDW phase develops at a higher temperature. 

We solve Eqs. (17) and (18) numerically, and show 
our results in Figs. m In agreement with an earlier 
work, 41 we find that for a narrow region in To values, 
there is a region in the phase diagram where SDW and 
superconductivity coexist. Specifically, superconductiv¬ 
ity emerges when To reaches some value denoted by r^. 
With further increase of intraband scattering, the SDW 
order is fully suppressed at some value T 0 = Fg Sdw \ As 
the ratio r^/To increases, both F^ and Fg Sdw ^ decrease. 
In Fig. [2] we plot the width of the coexistence region 
(rg Sdw) — rg SC ^)/27rT c0 at T = 0. Thus, we conclude that 
the coexistence region remains quite robust with respect 
to the interband scattering, and it only vanishes when 
both scattering rates become comparable, Tq ~ 



FIG. 3: (Color online) Phase diagrams in the (Fo,T) plane 
for T s o = 1.7T c o. In panel (a) we show the variation of the 
critical temperatures with Fo for = O.lFo, while panel (b) 
shows critical temperatures for F^- = 0.3Fo. The insets show 
the variation of the pairing amplitude and magnetization with 
disorder, evaluated at T = 0.01T c o. We use these data to 
evaluate the variation of the London penetration depth with 
disorder and temperature. 
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III. LONDON PENETRATION DEPTH 


In this section we solve the Eilenberger equation (13) 
and use the resulting correction to the Eilenberger func¬ 
tion to compute the penetration depth as a function of 
disorder and temperature in the lowest order in A. 


A. Solution of the Eilenberger equation in an 
external magnetic field 


The field-induced correction to the Eilenberger func¬ 
tion (14) is given by the solution of the following matrix 
equation 



S w T3P3(Jo, & 1} (R>n) 


ev F n ■ Af 0 /3 3 (Jo,^ 0) 


(19) 


which is found from (13) by keeping terms linear in the 
vector potential. The function C/), 1 ■* (R, n) must also sat¬ 
isfy the following condition, which results from the nor¬ 
malization of the full Eilenberger function (fl4|): 


gW . g(J) (R, n) + gW (R, n) • ^ 0) = 0. (20) 


We look for a solution of this equation in the following 
form 


^ 1} (R, n) = g£\ R, n)f 0/ S 3 <Jo - R, n)f 3/ 5i<T 2 

- ^(R, n). 


( 21 ) 


The matrix form for the first two terms follows from solv¬ 
ing Eq. (191; first in the limit when A = M = 0, and then 
for M = 0. In order to find the matrix structure of the 
third term we use condition (20), which can only be ful¬ 
filled for 


4 1} (R> n ) = S&T2P20- 1. 


( 22 ) 


Using condition (20) we obtain for the functions g£\ 
fu \ and sL 1} : 


with 


(R-i n ) = ev F n ■ A(R), 

Z(jJ 

fvH R,n) = -?^-ev F n- A(R). 

Z(jJ 

s^(R, n) = — ev F n ■ A(R), 

Zijj 


f 2 

J U) 


Zu J — (“b fh t9ui)9u) “b (A fr sfu)fu) 

+ ( M + 


(23) 


(24) 


where T s = T 0 — T^. Equations (23) and (24) consti¬ 
tute the perturbative solution of the Eilenberger equation 

pl. 


FIG. 4: (Color online) Inverse square of the London pen¬ 
etration depth (dimensionless units) as a function of intra¬ 


band scattering rate To evaluated at T = O.OlTco. In the co¬ 
existence region A)) 2 grows linearly with A in contrast with 
the clean case where A£ 2 ~ A 2 , Eq. (311. 




FIG. 5: (Color online) a) Temperature dependence of A£ 2 
for various values of To. b) Temperature dependence for 
5Xl(T) = Xl(T) — Al( 0) at low temperatures. The inset 
shows the same data as in the main panel-(b) but plotted 
in log-log scale. The interband scattering rate is fixed at 

= o.ir 0 . 
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B. Effect of disorder on the penetration depth 


In the clean limit, T t = 0, it is easy to show that 


The expression for the current density in terms of the 
original Green’s function is given 


eT \ ^ 

3 E 


e? 2 p 


2m ' J (27r) 2 

UJ n V 7 

Ne 2 A 


pTr 


fspoo-oG^lR.p) 


(25) 


where the last term guarantees the gauge invariance of 
the normal state. We can now use the same approxima¬ 
tion that we have already employed in our derivation of 
the Eilenberger equation: since the main contribution to 
the current comes from a narrow energy region around 
the Fermi surface, in the integral above we approximate 
p k pj?n. It follows that 


3 = 


—ie'KvppT ■ 
8 m 


, ro ( 03 do^ 1 ) (R, n) 


(26) 


where we used the definition ( 121 ; angular brackets de¬ 


note averaging over directions of unit vector n. Using 
our solution of the Eilenberger equation together with 
the equations for the current we obtain at an inter¬ 
mediate step 


j = -Q A, Q = -- 


2 c 




(27) 


where the auxiliary function is defined by Eq. (24). 


This expression can be significantly simplified by using 
the mean-field equations together with the normal¬ 
ization condition. Indeed, one observes that first and 
third terms in can be combined as follows 

(cu n T r tQLo^dto i(Al T — 

(kfo T ^ t9uj)9ix (ifo/pk, fos^poj) 

9u> 

s 2 H f 2 

{ui n + Ttg^g^ — (uj n + T t g u )— = (w„ + T t g u ) -—, 

9u 9ui 

(28) 

and then (w n +r t ^)(l + / 2 )-?;(A-ir s / w )/Z^ = (w n + 
Ttgui)- Consequently Q can be brought to the form 


„2„,2 


Q = 


2 c 


lT E 


f 2 

J cudco 


r t gu. 


(29) 


Curiously, the function does not enter explicitly into 
the final expression for the current. Thus, for the London 
penetration depth we have 


AZ \T) = 


2„,2 


ve v 


:2 ^E 


fu>9x 


u n + r tgui ’ 


(30) 


which is the main result of this paper. Next we analyze 
various limiting cases. 


AZ =AZ 2 (0) 


A 2 


M 2 + A 2 ’ 


(31) 


in agreement with earlier studies. ^ —I Let us now analyze 
the Matsubara sum in p0| ) in the limit of low temper¬ 
atures and for weak interband disorder, T, < Tq. In 
that limit, using Eqns. for the function f we find 
foj « iAg^/oo. Next, we set T —>■ 0, and convert the 
frequency summation into an integral over the variable 
x = uj/ A. The resulting expression for the Matsubara 
sum in (30) has the following form 


2 ^E 


x[x + "fg(x)] 1 dx 


1 


M 2 x 2 

A 2 [x +27 g{x)] 2 


3/2 ’ 


(32) 


where we have introduced the parameter 7 = fo/A. For 
small enough values of A, such that 7 /g> 1, it follows 
that for the moderate range of x ~ 0 ( 1 ) we can sim¬ 
plify x + 2y g{x) ss 2"/g(x). Furthermore, the dominant 
contribution to the integral ( |32[ ) comes from the region 
of x where g{x) ~ x. Thus, for sufficiently small A, we 
approximately obtain for the integral 


AZ 2 * AZ 2 (0) 


A 




2 + 4T 2 


(33) 


In the opposite limit, and still at zero temperature, an¬ 
other analytical result for Al as a function of the scat¬ 
tering rates and A can be derived (see Appendix B for 
details) 


AZ 2 = A 


- 2 ( 0 ) 
7tt( 3 + 8y 2 ) 

12 7 I 


tt(7^ + 27^7 2 + 4 y 3 ) 
127 s 4 

_ (7tt + 47s) arccos( 27 s 

~ 4 7 s 2 


(34) 


where we introduced the following parameters for brevity 
7 S = 2Ts/A, and 7 ^ = 2T W /A. Note that when inter¬ 
band scattering becomes negligibly small, then the term 
of Eq. ( |3lj ) in square brackets is proportional to A/To, 
in agreement with our estimate, Eq. (33), taken in the 
same limit. 

Our complete numerical analysis of equation (|30|) con¬ 


firms our asymptotic analytical expressions, and in par¬ 
ticular the estimate (33). In fact, we find that this be¬ 
havior persists for much larger values of A ~ T t . Lastly, 
comparing this result with the corresponding expression 
for the London penetration depth in the clean limit (31), 


we conclude that disorder has a crucial effect on the de¬ 


pendence of \ L " on both A and M. 


In Fig. [ 4 ] we show the variation of \ L 2 with dis¬ 


order, computed using Eq. (30) together with A(T 0 ) 
and M(To), which in turn have been computed self- 
consistently and are shown in Fig. [3] Furthermore, in 
Fig. I we show the temperature dependence of \ L 2 (T) 
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FIG. 6: (Color online) Plot of A)) 2 on A across the phase dia¬ 
gram at T = O.OlTco. Remarkably, the position of the line for 
the purely superconducting state depends on the strength of 
the interband scattering rate. This is expected as the super¬ 
conducting order parameter is strongly suppressed for higher 
values of F„/Po. 

for various values of To across the phase diagram [see 
inset (a) in Fig. [3]. Finally, in Fig. [b] we show the de¬ 
pendence of A on A^ 2 evaluated at T = 0.01T cO . One 
immediately observes that combining measurements of 
the London penetration depth in the coexistence region 
with those made in the superconducting state should, in 
principle, allow one to obtain an estimate for the ratio 
between the interband and intraband scattering rates. 
Indeed, for moderate values of r^/To, superconductiv¬ 
ity is strongly suppressed, leading to lower values of A 
in the superconducting state compared to those in the 
coexistence state. 


IV. DISCUSSIONS AND PERSPECTIVES 

In this paper, we obtained the phase diagram of doped 
iron-pnictide superconductors and calculated the mag¬ 
netic penetration depth at different temperatures under 
the assumption that doping introduces disorder but does 
not affect the band structure. In several limiting cases 
we have been able to reproduce previously known results. 

Our main finding concerns the behavior of Xl in the 
coexistence phase that has not been systematically ana¬ 
lyzed before in the presence of disorder scattering. Our 
modeling shows that starting from the overdoped side, 
A£ 2 grows with the reduction of the scattering rate in¬ 
duced by doping up to an optimal doping where super¬ 
conducting order is maximal. Once the system enters the 
coexistence phase, the inverse square of the penetration 
depth, which is proportional to the superfluid density, ex¬ 
hibits a kink followed by a sharp falloff. In that region, 
disorder primarily affects magnetic order rather than su¬ 
perconductivity. In sharp contrast with the clean case, 
where X'J 2 oc A 2 near the end-point of the superconduct¬ 
ing dome from the side of the pre-existing SDW, within 


the disorder model we find a completely different scaling 
law A£ 2 oc A. Another important observation concerns 
the low temperature dependence of 5Xl(T). The log-log 
plot presented in Fig. (|5]o) suggests either a power law 
behavior of SXl on temperature, SXl oc T a , with rather 
high power exponent a > 5, or exponential dependence 
SXl oc e~ w / T . The exponential dependence indicates the 
presence of the gap W in the spectrum of electron states 
near the Fermi surface for parameters of curves presented 
in Fig. ([5]). It is plausible, however, that for stronger 
electron spectrum becomes gapless at low temperature 
and SXl{T) exhibits a power law with a ~ 2P^ 

Several extensions of the presented model are in order 
to improve the com pariso n with experiments. First, one 
could treat the bancP^EIl and disordei®^ models on an 
equal footing in order to study the observed anisotropy in 
A l- Second, one could consider the extension of the pre¬ 
sented formalism beyond the Born approximation, which 
might be necessary for an accurate interpretation of the 
low-temperature data. Third, one could account for dif¬ 
fusive scattering from the surface of the superconductor. 
That would complicate the calculation of the penetra¬ 
tion length considerably, since one would have to work 
with an integral Milne equation, instead of a differential 
London equation, which governs the distribution of the 
magnetic field in superconductors and consequently de¬ 
termines the precise value of Xl in the nonlocal limit 
Fourth, our analysis, so far, has been restricted to the 
mean-held level. A most intriguing recent experimen¬ 
tal observatior l 48 l 49 l is an apparent sharp peak in A^ ob¬ 
served in isovalently P-doped BaFe 2 As 2 at nearly zero 
temperature around the optimal doping. This effect was 
attributed to quantum critical fluctuations of the SDW 
order at the onset of the transition into the coexistence 
phase I52H52] jt j s 0 f apparent theoretical and experimen¬ 
tal necessity to investigate to what degree such quantum 
effects are robust against disorder scattering. On a tech¬ 
nical level, this would require the inclusion of magneti¬ 
zation fluctuations into the existing formalism. We note 
that such a generalization has already been carried out 
in the context of thermal magnetic fluctuations, which 
are relevant for the interpretation of specific heat dataP^ 
Finally, when analyzing the quantum critical behavior of 
the superfluid density in FeSC compounds, it might be 
useful to use the results obtained in the context of cuprate 
superconductors.^ In the vicinity of the quantum criti¬ 
cal point (QCP), generic scaling analysis indicates that 
the superconducting critical temperature should vanish 
as T c oc 5 ZU , where 5 = \x — x c \ measures the devi¬ 
ation in doping from the QCP, while z and v are the 
quantum dynamical and correlation length exponents. 
At the same time, the superfluid density should scale 
as n s oc X 1 [ 2 oc S < ' z+d ~ 2 ' >v , where d is the dimensionality 
of the system. When combined together, these two scal¬ 
ing laws predict that there should exist a precise relation 
dlnT c /dlnn s = z/(z + d — 2). In the two-dimensional 
case there should exist a linear relation between T c and 
n s . For FeSCs, each of the end-points of the coexistence 











phase represent a QCP and, consequently, establishing 
the relation between T c and n s will provide new informa¬ 
tion about superconductivity in these complex materials. 
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Appendix A: Derivation of the Eilenberger equation 

We start by writing down the equations of motion for 
the matrix Green’s function using an imaginary time rep¬ 
resentation, while ignoring the disorder potential for the 
time being: 

- ^-G(riTi;r 2 T 2 ) - H{yi) ■ G(y 1 t]y 2 t 2 ) 

OTl 

= I6(t 1 -t 2 )5{ r!-r 2 ), (Al) 

d - 

—G(riTi;r 2 r 2 ) - G(y 1 t 1 \y 2 t 2 ) ■ H( r 2 ) 
ot 2 

= I5(ti - r 2 )(5(ri - r 2 ). (A2) 

Next, we use mixed space-momentum, keeping the 
center-of-mass coordinate R = (ri + r 2 )/2, but making 
the Fourier transformation with respect to the relative 
coordinate r = iq — r 2 . In the Matsubara frequency rep¬ 
resentation, we find: 


i^n^oPO^O + £p — 


ip • V 


R 


T3P3^oG w (R,p) 


p 2m 

+v • A(R)f 3/ 5 0 (JoG tJ (R, p) - R m f(R)G aJ (R,p) = I(A3) 
iUnG^i R, p) + Gu)(R, p)t3/3 3 <5-o + 

G U (R, p)f 3/ 5 0 (7ov • A(R) - G W (R, p)R mf (R) = /, (A4) 


where we redefined the vector potential -A —> A for 
brevity, and we used Eqns. ® and Q. We write 
^ ps vf n, v ps ^ 11 , n = 27 and using Eq. (Il2|) we 
mult iply Eq. ( |A3| by t 3 p 3 ctq from the left, and multiply 
Eq. (A4| by the same matrix from the right. Then, we 


subtract the second equation from the first. Lastly, we 
use T 3 / 931 J 0 • ' 7 p 3 d'o = I where necessary, and integrate 
both parts over the absolute value of the momentum, 
which allows us to use the Eilenberger Green’s function 


(121. We thus find 


iiVnGui R- n ) “ iw n f 3 p 3 CTo^(R, n)f 3 p 3 a 0 
+u F n ■ (-iV R )f 3 p 3 (To^(R, n) 

+v F n ■ A(R) f 3/ 5 0 (Jo^(R, n) - f 3/ 5 3 (Jo^(R, n)f 0y 6 3 CTo 

+T 3 p 3 aaQ U] (R., n)i7 mf (R)f3/5 3 clo - 
-f 3 p 3 a 0 H mi (R)f 3 p 3 a a g u (Il, n) = 0. (A5) 

This equation can be written in a compact form if we 


multiply it from the left by f 3 p 3 ao. Then, Eq. (13) from 
the main text follows, where we have included the effects 
of disorder by trivially writing the self-energy correction 
to the mean-field Hamiltonian. In terms of the Eilen¬ 
berger Green’s function it is given by 

C d(j) 

E^(R) = -ir 0 f 0 p 3 a 0 T 3 p 3 a 0 / -^£ W (R, n)f 0/ 03<T 0 

0 J d ^Gu(R;n)T 1 p 3 ao. (A6) 


—*17 fi p 3 (Jof 3 p 3 (Jo 


Appendix B: London penetration depth for M <C A 


In this section we will derive an expression for the Lon¬ 
don penetration depth at low temperatures, assuming 
that the SDW order parameter is much smaller than the 
superconducting order parameter, M< A. We start by 
writing Eqns. © where we replace / —> if and s -A is: 

Ag w = (wn + 2T w )f u , Mg u = [uj n + 2r t )s„. (Bl) 


The functions in (Bl) satisfy t he n ormalization condition 
9u> + ft + s ui = 1- From Eqs. (Bl) it follows that 


A- M/7 = -2 T 0 f u s u 


(B2) 


We can now eliminate from this equation by using 
the normalization condition, which yields the following 
equation for g u : 


1 — (7 2 — f 2 

JLU J UJ 


1 - 


A 

w. 


V 1 - 9l - fl 


1 2 


4r 2 

= m (B3) 


Solving this equation for g^ one obtains 




M 2 f* 


(2 T 0 / w + A)2 ■ 

Next, we consider the following integral 


Q(A,M) = - / 
Jo 


0 w + r t gu) 


(B4) 


(B5) 


which in a way determines the penetration depth. The 
idea is to replace the integration over ui with an integral 
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over fa. To do that, we employ Eqs. (Bl|, (B2), and 
(B41. It follows then that 


Q(A,M) = [ 
Jo 


fudf u 


[A 9l 


I o 9u{ A + r s /tj) 
+(A-2T w / a Wl + 


(A + 2T 0 fu) 3 J\ ’ 


(B6) 


and g u is a functional of f u , Eq. ( |B4[ ). Clearly, for T t = 0 
we find Q(A,M) = A 2 /{M 2 + A 7 ). When M = 0, the 
expression for Q{ A, M = 0) simplifies to: 


Q(A,0)= f 

Jo 


(A - 2r 7r / 3 )/d/ 


0 yfof(A + 2r s /)' 


(B7) 


This integral can be evaluated exactly, and it gives 
Eq. (34) from the main text. Lastly, one can also ex¬ 


pand Q( A, M) in powers of M/A to derive the correction 
to the penetration depth due to the development of the 
SDW order in the superconducting state. The resulting 
expression, however, is too cumbersome to show here. 
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